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Abstract 

Motivated by recent findings due to Wiegmann and Zabrodin, Faddeev 
and Kashaev concerning the appearence of the quantum U q {sl{2)) symmetry 
in the problem of a Bloch electron on a two-dimensional magnetic lattice, we 
introduce a modification of the tight binding Azbel-Hofstadter Hamiltonian 
that is a specific spin— S Euler top and can be considered as its "classical" 
analog. The eigenvalue problem for the proposed model, in the coherent 
state representation, is described by the S— gap Lame equation and, thus, is 
completely solvable. We observe a striking similarity between the shapes of 
the spectra of the two models for various values of the spin S. 
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1 Introduction 



The quantum mechanics of free electrons on two dimensional lattices, in the presence 
of a homogeneous and transverse magnetic field, ("magnetic lattices") leads to the 
discovery of a host of physically and mathematically fascinating problems. 

This subject has a long history, starting with the pioneering papers of Harper, 
Azbel, Zak and Chambers, Hofstadter and Wannier 0, 0, § g, §, §. With the 
discovery of the quantum Hall effects |7], |8|, ||, a large number of very interesting 
theoretical papers appeared, which deal with the quantum mechanical explanation 
of the Hall conductivity plateaus [LI] . 



More recently, in the work of Wiegmann and Zabrodin []n|, it was found that 
the eigenvalue problem for the the Hamiltonian of the one electron lattice problem 
(for rational magnetic flux per plaquette), henceforth called the Azbel-Hofstadter 
(AH) Hamiltonian, can be written as a q— difference quadratic equation, using the 
quantum group U q (sl{2)), known as the "Jimbo" deformation of SU(2), which ap- 
pears naturally in this context. Their analysis leads to the algebraic Bethe Ansatz 
equations for the roots of its polynomial solutions, at a particular point of the Bril- 
louin zone of the square lattice. Subsequently, Faddeev and Kashaev proved that 
this symmetry exists at all points of the Brillouin zone for square and triangular 
(anisotropic, in general) lattices and they provided the corresponding Bethe Ansatz 
equations ]Rj] . 

However, it has not yet been possible to solve the Bethe Ansatz equations thus 
obtained. Direct numerical solutions were found, in certain cases by Hatsugai, 



Kohmoto and Wu [jig] ; but the necessity of a systematic approximation scheme 
remains an open issue. 

In this paper we propose a modification of the AH Hamiltonian through a specific 
spin— S Euler top, which has the merit of being completely solvable. Indeed, in the 
coherent state representation, the eigenvalue problem for this Euler top is described 
by an S— gap Lame equation. Numerical comparison of the spectra of the two models 
reveals a striking similarity for their shapes. 

In section we establish a uniform notation, recalling, at the same time, the 



salient results of ref. [12, 151. In section RL we write the AH Hamiltonian in terms 



of the generators of the Cartesian q— deformation JT9[ of SU(2) and set the stage 
for a model Hamiltonian, 7ii, which is that of an "Euler top" under the ordinary 
SU (2) group. In the process we will establish connections between the two different 
q— deformations of SU(2), a problem that is interesting in its own right. 

In section [| we discuss the eigenvalue problem for the AH Hamiltonian and 
present the recursion relations for the eigenvectors and the eigenvalue equation in a 
compact 2x2 matrix form. The case E = admits an explicit solution [L6|] . 

In section [| we explore the symmetries of the classical Hamiltonian, Ti,±. We 
provide explicit recursion relations for the components of the eigenvectors and for 
the eigenvalues. We also show how it is related to the S— gap Lame equation, using 
the coherent states of SU{2). 

In section |6] we provide numerical evidence that the spectrum of the AH Hamil- 
tonian may indeed be meaningfully approximated by that of Tii . 
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We end with our conclusions and a discussion of directions of further inquiry. 



2 The Quantum Group sl q {2) and the AH Hamil- 
tonian 

The Azbel-Hofstadter (AH) Hamiltonian is a tight-binding model for a single Bloch 
electron on a two dimensional lattice and in the presence of an external, homogeneous 
and transverse magnetic field P| [3| |5|, |5|, 
The AH Hamiltonian is 

^ah= t nrn e iAnm clc m (1) 

<m,n> 

where t nm are the hopping amplitudes, , c n creation and annihilation operators 
for the electron at site n = (n x , n y ) G Z x Z and A nm the line element, viz. 



e 

C J n 



in 



A-dx = A nm (2) 

such that, through each elementary plaquette, we have a flux $, 

Y[ e [Anm = e^^° (3) 

plaquette 

where <fi = $/$o = 2ttM/N with M and A coprime integers and <3>o the flux 
quantum. 

For such values of the magnetic flux the Hilbert space becomes an infinite number 
of identical A— dimensional copies, due to the existence of an infinite number of 
magnetic translations, that commute with Hamiltonian 7Yah- In the case of a square 
lattice the dynamics is restricted to a Ax A square, magnetic, lattice. In the Landau 
gauge, 

Ac = A nx>nx+l = A y = A ny>ny+1 = $ • n x (4) 

the Bloch wavefunction ip(n x ,n y ) = exp(i(k x n x + k y n y ))ipk{n x ) , where (k x ,k y ) G 
[0, 2tt) x [0, 2vr) and 

M^x + N) = xp k (n x ) (5) 

satisfies Harper's equation 

t x (e ih ^ n+1 + e^n-i) + t y (e^y +n ^ + e" 1 ^^)^ = E^ n (6) 

where t x ,t y are the hopping parameters and ip n = ip^{n). 

We shall restrict our discussion to the case of an isotropic square lattice, t x = 

ty = l- 

Harper's equation can be written in matrix form 

HV = (7) 
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where \I/ = (ipi, . . . , iPn) t and 

H = e ikx P + e-^p- 1 + e lky Q + e'^Q- 1 (8) 
The matrices Q and P are 

(9) 

P M = 4-1,2, k, I = 1, . . . , N 

with uo = exp(27riM/iV) and all operations are performed mod N. 

The matrices P and Q generate the magnetic translations and the finite Heisen- 
berg group through the Weyl commutation relation 

QP = uPQ (10) 

The Heisenberg group elements 

J r<8 = u r - s/2 P r Q s (11) 
provide a projective representation of the two dimensional translation group mod 

\ m q 

Jr,sJr',s'= ljji ' r S T Jr+r',s+s' (12) 

The factor 1/2 in the exponents of the previous relations is defined as s + 1 for 
N = 2s + 1 and s an integer, while, for N even, it is defined as 0.5 . The matrices 
J r ^ s are traceless, unitary and they have period N, viz. 

J r s J—r, — s 

(13) 

TN — 1 
<Jr,s ~ X N 

The matrices Q and P are related through the finite Fourier transform matrix F, 

F k>l = -L^^-DC-i), k,l = l,...,N (14) 



by 

F-P = Q-F (15) 

Because of the symmtery 

Jl. ■ P ■ Jr,s = LO- S P 



JL ■ Q ■ Jr,s = 



(16) 



the Brillouin torus [0, 2n) x [0, 27r) is reduced to a smaller one [0, 2%1/N) x [0, 2irl/N), where 
I = M mod N (recall that the flux is 2tcM/N). 

In ref. jl2| Wiegmann and Zabrodin made the important observation that, for 
the midpoint of the Brillouin zone, 

cos Nk x + cos Nk y = (17) 
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the spectrum of Hah is determined by the roots z m of polynomials P(z) of degree 
N — 1, that interpolate the wavefunction 

^ n = P(q%), n=0,...,N-l 

(18) 

q = e ™ M ' N 

The energy spectrum is given by 

N-l 

E = iqS(q - q, 1 ) £ z m (19) 

m=l 

and the roots z m satisfy the algebraic Bethe Ansatz equations 

2 

lO ~ ^fc _ -Q % Z k — Zm ^Q) 

^O^fc — 1 mf L k Zk — %Z m 

where k — 1, . . . , N — 1 . 

Behind these findings is the quantum group sZ go (2), which is a deformation of 
the Lie algebra sl(2) jl7| 

2J ^-2J (21) 

% - % 

where J±, Jo are expressed in terms of the magnetic translations J r ^ s and the Hamil- 
tonian, in a specific, "quasi-Landau" gauge, is expressed, in terms of the J± as 

H A n = i(qo-q 1 )(J-±J+) (22) 

at the midpoint (s) of the Brillouin zone. 

In ref. |15| Faddeev and Kashaev observed that the AH Hamiltonian is equivalent 
to the three site chiral Potts model whose quantum group symmetry is known 
to lead to the algebraic Bethe Ansatz equations over specific Riemann surfaces. Thus 
they were able to generalize the result of Wiegmann and Zabrodin to arbitrary points 
of the Brillouin zone, as well as to anisotropic square and triangular lattices. 

3 The Cartesian deformation of sl(2) and the AH 
Hamiltonian 

Some years ago a new deformation of the sl(2) algebra, in the Cartesian basis, was 



proposed, which has a cyclic symmetry for the generators |19 



qJiJ 2 


- q~ 


J2J1 — 


(q 2 


-t 2 : 






- q~ 


-V3J2 = 


(q 2 
(q 2 


-q~ 2 : 


Jl 


qJsJi 


- q 




-q- 2 : 





(23) 
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The representation theory for real values of the deformation q was studied by Fair- 
lie [1J| and, in more detail, by Zhedanov et al, who pointed out that this algebra 



generates the properties of a class of Askey- Wilson polynomials. Recently in ref. |2C 
the representations for q a primitive root of unity of order N were constructed and 
classified. It is known that, in this case, the standard ("Jimbo") deformation of 
sl(2) has irreducible representations of all dimensions smaller or equal to the order 
of the root. These representations depend on three complex parameters and are of 
two types: type A irreps have a classical si (2) analog, while type B are cyclic (i.e. 
not ladder) of order N pi 



The Cartesian algebra, sl^ (2) has one Casimir element 

C 2 = (q + q- 1 )^ + Jl) - { J 3 , J 3 } (24) 

where (q 2 — q~ 2 )J 3 = q~ x J\J 2 — q^Ji- In ref. [^2], ^3|, [25|], using results from Finite 
Quantum Mechanics, it was realized that the AH Hamiltonian,at any point of the 
Brillouin zone, can be written as the anticommutator of the two operators J\ and 
J 2 , 

Hah = —^Uu J2] (25) 
q + q^ 1 

where J±,J 2 and J3 realize an N— dimensional representation of slg (2) and may be 
expressed in terms of the generators of the Heisenberg group, J rjS through 

L = p ia J -\- p-' la T 

°1 °m,—m ~ u —m,m 

Jl = &^Jm,m~\~& ^J—m,—m (26) 
J 3 = e ikyQ + e ~ik y Q-l 

where q = u m \ u = e 2ni ^ N ,m = (l/2)mod N = (N + l)/2 and a = {k x - k y )/2, 
p= (k x + k y )/2. 

The unitary representation of sl^(2) constructed above is irreducible with Casimir 
element 

2ttM . , 

C 2 = 8 cos -j^™ (27) 

and 

J 3 = e ik x p + e -ik x p -i ^ 

The generators Ji, J 2 , J3 can be cyclically permuted at the Brillouin point (k x , k y ) = 
(0,0) by three N x N unitary matrices, Ui, U 2 , U 3 

U-\] 2 V X = J 3 

U 2 l J 3 U 2 = J± (29) 
U 3 l J x U z = J 2 

Theses matrices leave the corresponding Jj's invariant. This situation brings to mind 
the rotation group in three dimensions, where n/ 2— rotations around the coordinate 
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axes cyclically permute the generators of the group. In the appendix [A] we construct 
explicitly the group generated by U\, U2, U3. 

Before concluding this section, we establish the relation of the iV— dimensional 
unitary matrices Ji, J2, J3 with matrices J± and Jo, which satisfy the (Jimbo) 
deformation of si (2), i.e. sl q (2). 

It is straightforward to see that, defining J± and J by 



J± i 1 (i7±m,±m Jtm.Tr 

q-q 1 



,m 2 



(30) 



Q = to mJ ° 

the Jimbo deformed algebra is satisfied. Then, introducing matrices 

T± = P-\J ± P (31) 
we establish the relation between J 1; J 2 , J3 and J± and Jo as 

Ji = ^^(g 2 J + + g- 2 J_-T + -T_ 
q + q 1 v 

</ 2 = -— ^ (g- 2 J + + g 2 J_ - T + - TJ ( 32 ) 
q + g v 

J 3 = g 2Jo +<T 2Jo 



where q = uo r 

Finally, we note that the constructed representation for J ± and J is of the cyclic 
type since we can check that eq. ( ^01) implies 

J^ = 2 x ljv (33) 



4 The eigenvalue problem for the AH Hamilto- 
nian 

Although in the literature the eigenvalue problem has been discussed in various 
contexts, we present here a compact method for the determination of the eigenvalues 
and eigenvectors, using the tridiagonal form of the AH Hamiltonian. This method 
is especially suited for fast numerical calculations for large N. 
The eigenvalue problem, in components, is the following 

e ifc ^fc-i + 2 cos Uk - 1) ^ + k^j ij k + e~ ik ^ k+1 = E^ k (34) 

where k — 1, . . . , N and we use periodic boundary conditions 

ip N+ i = ipi (35) 
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for all I. We define homogeneous variables Zf. 



z k = — — (36) 



and we use the "Mobius" notation 

. I ■ 7 — 

cz + d 

to rewrite the equations fl34|) in the form 




(37) 



E-2cos((k-l)^ + k x ) -e ik y 



zk = ^ .^^,„,^ . y zki (38) 

for k = 1, . . . , N. The last equation, k = N, gives, by iteration, the relation 

^g(*-M(*-/>f+**) -f )}.,„ (39) 

Since z^ = z , we deduce the characteristic equation, viz. 



det 




£-2cos((£;- l)^ + k x 









| — 12x2 





(40) 



and we may then compute the components of the corresponding eigenvector, by 
recursion, starting from z^, which is a fixed point of the Mobius transformation, 
eq. ©. 

It is known @] that the characteristic polynomial takes the form 

P N (E) + (-) N x 4 (cos Nk x + cos Nk y ) (41) 

where Pn{E) is a polynomial of degree N in E and coefficients that do not depend 
on (k x ,k y ). The study of the structure of the gaps between the eigenvalues is a 
very interesting problem in functional analysis and noncommutative geometry (cf. 
ref for a recent review) 

It should be noted that the existence of a zeromode for eq. ([E]) depends on 
(k x , k y ). We observe that, for iV even and r = AT/ 2, we obtain 

jln A RJr,r + n A R = (42) 

for any point of the Brillouin zone. This implies the existence of a reflection sym- 
metry for the spectrum, i.e. for any eigenvalue E there exists an eigenvalue —E. 
This does not necessarily imply the existence of the eigenvalue E — 0; if it exists, 
however, it must necessarily be doubly degenerate. 

For iV odd the reflection symmetry, E <->• —E is approximate. 

On the other hand, for k x + k y = tt, the Hamiltonian anticommutes with the 
discrete Fourier transform, which implies that, for any N, the reflection symmetry 
is realized and, for iV odd, these points belong to the midband, eq. ( |TTD and the 
zeromode always exists. 

In the context of the quantum group U q (sl(2)), the eigenvector corresponding 
to E = is a g—Askey- Wilson polynomial |L2| and the structure of its roots as a 



function of the magnetic flux has been studied by Kohmoto et al. |jl6 
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5 The symmetries of an Euler top 

We introduce now the classical analog Hamiltonian 

H 1 = j i {S 1 ,S 2 } = l i (Sl-S 2 ) (43) 

in the spin-S* representation, with Q = 2S + 1, where Si are the standard SU(2) 
generators. 

We propose to study the eigenvalue problem of Tii with the hope to gain some 
intuition, which may be useful for the real problem (i.e. that of the AH Hamiltonian). 
The matrix elements of Tii are 



4i 

where 



(WiW = ^ (aA',fc+ 2 - a A,fc- 2 ) ( 44 ) 

afc = a k s-iak-s (45) 

and 

a m = y/S(S + l)-m(m + l) (46) 

with k, k' = 1, . . . , 2S* + 1. The structure of the matrix shows immediately that the 
even-numbered components decouple from the odd-numbered ones for any value of 
S. 

The three operators 1Zi >2 = e~ 17r,Sl ' 2 and 1Z 3 = e 1 ? 53 . will prove useful for the 
decomposition of the 2S + 1 -dimensional eigenspace in convenient subspaces. In 
components they read as follows 

i n i)k,k' = emS h+k',2S+2 
)k,h 

They satisfy the following commutation relations 



(Ks) k , k > = ^ 8 (-i) h -%y (47) 



n{n 2 = {-iy s n 2 K 1 = (K 3 y 

n 3 n x = n 2 n 3 

(K 1:2 ) 2 = (n 3 r = (-l) 2S (48) 

From the above it follows that ^.1,2,3 anticommute with the Hamiltonian and we can 
construct two projectors, that commute with the Hamiltonian, 

V ± = l(l± e^ s (TZ 3 ) 2 ) (49) 



and 



Q± = l(i±e-^ s n 3 n 2 ) (50) 
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V± project on the subspaces of the odd- and even-indexed components of the eigen- 
vectors, while Q± on the positive and negative energy eigenspaces. The Hamiltonian 
TCi thus may be written as a direct sum of Hamiltonians 7i± 



n = n+®H 



(51) 



where 



u± = nv. 



(52) 



and the dimensions of the corresponding Hilbert spaces are, for integer S, S + 1 

and S respectively while, for half-integer S, they both have dimension S + — . In 
components 



It is straightforward to numerically diagonalize the Hamiltonians 7i± and compute 
the eigenvalues and eigenfunctions. The anticommutation relations ( fl8|) imply a 
spectrum antisymmetric about E = for any value of S and for both subspaces. 
The general structure of the eigenvectors is qualitatively similar to that of the one- 
dimensional harmonic oscillator and their phase structure (real/imaginary compo- 
nents appear symmetrically or antisymmetrically up to factors of ±1 or ±i) may 
be deduced from the antidiagonal structure of the projectors Q±. Analytically, it is 
possible to completely describe the E = case and provide recursion relations for 
the E 7^ cases. 

We introduce the standard notation for the Mobius transformation 



It is easy to show that the eigenvalue problem 7i\ty) = E\^f) takes the following 
form in component notation (where £ = 2iE) 



a2/c+A',fc+i ~ a 2fc-A',fc-i 



(53) 




(54) 



a 2fc-1^2fc-l + a 2fc+lV'2fc+3 
— 32fc'^2fc + 32fc+2'?/'2fc+4 



(55) 



corresponding to the odd and even subspaces. Defining 



1p2k-l 
4>2k+2 
^2k 



(56) 



we solve eqs. (|55|) by the Mobius transformation 




(57) 
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where k = 1, . . . , S — 1 and z\ = £/a\. The last equation in this subspace takes the 
form 

°=(t *V)- zs (58) 

which is the characteristic equation for E in the odd-numbered sector. Similarly, 
for the even components one has (w\ — £/az) 

{n(al a r)}- (») 

and the corresponding characteristic equation 

= ( t 325 - 2 ) ■ -5-! (60) 

The above analysis becomes explicit for E = 0. There are several possibilities, 
depending on the values of S. 

• integer S: 

E = is an eigenvalue only for S = 4k + 1, 4k + 3 (non-degenerate) and 
5 = 4 A; + 2 (doubly degenerate). 

The explicit form of the corrresponding eigenvector \ty b E=0 ) = (^i, . . . ,ip2S+i) 
is 

, ll m =0 a 4m+2 , / c1 \ 

^4n+2 = -j^j 1 ^2 (61) 

llm=l a 4m 

for the non-degenerate case and identical to that of the half-integer S (see 
below) for the doubly degenerate case; all other components are zero. 

• half-integer S: 

There are two possible cases to consider: S = (4k + l)/2 and S = (4k + 3)/2. 
In the first, E = is an eigenvalue, that is doubly degenerate, one belonging 
to the (+), the other to the (— ) subspace. 

The two eigenvectors, corresponding to E = are l 1 !^ 1 )) = (i/j^, . . . , V4s+i) 
and |*( 2 )) = (#,... 

1>. 



,(!) - 


rrn-1 
llm=0 


34m+l 


HP 


4n+l — 


Ilm=l 


34m- 1 


,(2) _ 


■prn-l 
llm=0 


34m+2 




4n+2 — 


n n 


1 34m 



^+2 = (62) 

llm=l a 4m 

For 5 = (4k + 3)/2 E — doesn't belong to the spectrum, as may be proved 
by explicit calculation. 

In the above relations the ipi,ip2 are normalization constants. To determine the com- 
ponents of the eigenvectors (zk,Wk) one has simply to evaluate the two-dimensional 
matrix products in eqs. (|57| , |59l) for each root £ of the characteristic equations (m,|0). 
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Another approach to the problem is provided by the coherent state representation 
of SU(2) where we write the generators as differential operators and the algebraic 
eigenvalue problem becomes a Schrodinger problem in a class of potentials. This 
connection has been used in the inverse way for classifying potentials that have 
a "quasi-integrable" spectrum, i.e. a finite number of eigenstates decouples from 
the rest and may be determined by finite matrix methods |23| . Even more 

recently, these same equations have appeared in the study of the correlation functions 
of Wess-Zumino-Witten models on the torus [F2^ . 

In fact TCi belongs to a family of Hamiltonians with the same spectrum 



n 



e iaS Wie- iaSs 



cos 2a • TCi + sin 2a ■ H2 



(63) 



where Hi = {Si,S 2 } and H2 = <Sf — <Sf. H2 is known to describe isotropic para- 



magnets in two dimensions and a similar Hamiltonian has been studied in ref. [£6 



Indeed, by a S2 rotation of angle 7r/2, one obtains the Zaslavskii-UTyanov Hamil- 
tonian H3 |26] in zero external magnetic field 



n. 



e -if<S 2e -ifS 3 



Hie 



i f- S 3 e if5 2 



^3 — <^2 



(64) 



In the coherent state basis, the SU{2) generators have the following form 

, d 



Si 

s 3 



S cos ( 



sin 1 



S sin + cos « 
d 



d 



(65) 



— 1- 



and the eigenvalue problem for H3 is 

'1 + cos 2 (j))-^- + (S - I) sin 2<t>4i + (E + S 2 sin 2 + S cos 2 
dm 2 dm 



$(0) = (66) 



The components ip m m = —S,...,Sof eigenvectors of H3, are related to the function 
$(0) by 



m= E 



Ipr, 



-s yJ(S-m)\(S + m)\ 
If we change variables, following ref. [[HJ, from (0, $(0)) to (x, ty(x)), defined by 



(67) 



y(x) = $(0(x))(l + cos 2 0(x))- s/2 

T = ^( 1 + cos2 ^)) 1/2 ^ 

dx a/2 d(p 



eq. (pq ), after a redefinition 



x = u- K(l/2) 



(68) 
(69) 
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where K(l/2) is the complete elliptic integral of the first kind |JIJ, becomes the 
S— gap Lame equation JHtJ 



d 2 E + S(S + 1) 
drf + 2 



5(5 + 1) 



sn 2 M 



V(u) = 



(70) 



where snu is the elliptic sine of modulus 1 / \[2 |JIJ . This equation has polynomial 
solutions in terms of elliptic sines and cosines and appears in the inverse scattering 
method for the KdV equation and the potential corresponds to the soliton solutions 
of the KdV. 



6 Numerical Results 

In this section we compute numerically the eigenvalues and eigenvectors of the AH 
Hamiltonian and the model Hamiltonian, TCi, for various dimensions. The size of the 
matrix representation of the AH Hamiltonian depends solely on the denominator of 
the magnetic flux, <fi = 2irM/N, where M and N are relatively prime integers. In the 
case of TCi, the size is = 25 + 1 and, as noted previously, there are two decoupled, 
invariant subspaces, of dimension 5 and 5 + 1 respectively. For simplicity, we shall 
compare the spectra, for N odd, of Hah and H\. 

Furthermore, our toy Hamiltonian does not depend on the Bloch momenta 
(k x , k y ) but, since the Brillouin zone has size equal to the flux, our approxima- 
tion should be better, for flux 2tt/N, at the point (k x , k y ) = (0,0), where the AH 
Hamiltonian has maximal symmetry (commutes with the finite Fourier transform 
matrix) and has both odd and even eigenf unctions, as is the case for 7i\. 

In fig. m we compare the spectra for the case N = 111 and in figure ||] we present 
the ground state eigenvectors for the two Hamiltonians. In comparing the spectra 
we have normalized them both to have the same numerical value for the ground state 
energy. In fig |3| we give the normalized difference of the two spectra for N = 111. 
In fig. H we display the maximum error for several values of N. 

7 Conclusions and Perspectives 

The main point of this paper is the striking similarity between the spectrum of 
the AH Hamiltonian and that of our Euler top model. For the moment this is 
an empirical observation, for which we still lack a physical appreciation from first 
principles. 

From the quantum group symmetry point of view the cyclic character of the 
U q (sl(2)) representations, which appear in this problem would forbid a direct classi- 
cal large A" limit. However, our numerical results indicate that our approximation, 
through the Euler top, gets better with increasing N. This probably suggests some 
kind of "analytic continuation" between the the Cartesian deformation of SU(2) 
and the "classical" SU(2) for the root of unity case. This is still an open problem, 
for which our results seem to indicate a promising direction of attack. 
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Figure 1: Comparison of the spectra of Hah (grey points) and Hi (black points), 
N — 111. 13 
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Figure 2: Comparison of the ground state eigenvectors of TC ah (grey points) and Hi 
(black points), N = 111. 14 
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Figure 3: Plot of E A n(n) - E H1 (n) x (Eau(1)/ E ni (l)). Note the striking difference 
between the values of the even and the orfd subspaces. 
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Figure 4: Log-Log plot for the maximum difference for N = 55 to iV = 145 in steps 
of 10 for the spectra E A n and E ni . Black dots are for the even subspaces and gray 
dots for the odd subspaces of Hi- 



Another problem is the dependence on the Brillouin zone parameters (k x , k y ) as 
well as on the flux = 2ttM/N; indeed we recall that our results hold only for 
M = 1 and it is an interesting question, for instance, to calculate the "Hofstadter 
butterfly" for the Euler top. 

We hope to return to these questions in future work. 
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A Quantum Rotations 

In order to construct the three matrices U\, U 2 and U3, mentioned in section [| we 
have to determine the inner automorphism group of the discrete Heisenberg group 
matrices J r ^ i.e. matrices U(A), such that 

U\A)J rtS U{A) = J {l , s)A (71) 

where r, s — 0, . . . , N — 1, for every matrix A 6 SL(2, Z^) |22], |2~3fl . It is obvious, 
from the definition of the generators J\, J 2 and J3, that we can construct matrices 
U\, U 2 and U3, if we find three 2x2 matrices Ai, A 2 and A 3 , which leave invariant 
the set of indices(m, — m), (m, m), (0,1) respectively. We check immediately that 
there exist three abelian subgroups of SL(2, Z N ), which do the job, generated by 

A 1 = 



A 



3 




We set Ui = U(Ai), with i = 1,2,3. Using the explicit forms for U(A), A E 
SL{2, Z N ) from ref. @, we find 



1 1 (73) 



(U 2 ) kjl = -^ojmd-W-u+V x _. 

(U 3 ) kjl = 5 kil co m ^ 2 
where k, I = 1, . . . , N and the symbol 

f 1 If 1N= lmod4 

\ -i J ~ \ -i, N= -lmod4 1 ; 

We note finally that the cyclic groups generated by Ui,U 2 and f/ 3 are of order N. 
The matrices Ui, U 2 and U 3 can be used to define specific discrete Askey- Wilson 
polynomials as the columns of the matrix 

V AW = U\U 2 (75) 

cf. the papers of Zhedanov and collaborators in ref. [[HJ. 
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